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Abstract 

The statistical properties of failure are studied in a fiber bundle 
model with thermal noise. We find that in agreement with recent ex- 
periments the macroscopic failure is produced by a thermal activation 
of microcracks. Most importantly the effective temperature of the 
system is amplified by the spatial disorder (heterogeneity) of the fiber 
bundle. 
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Material failure is a widely studied phenomenon not only for its very 
important technological applications but also for its fundamental statistical 
aspects, which are not yet very well understood. Many models have been 
proposed to give more insight into the statistical analysis of failure. Among 
the most studied ones we can mention the fuse and the non-linear spring 
networks [1], 0, f|, |5| which can reproduce several features of crack pre- 
cursors, experimentally observed in heterogeneous materials subjected to a 
quasi-statically increasing stress || 0, |8|, |9|, [10). Specifically the power law 



behaviour of the acoustic emission observed in several experiments close to 
the failure point. 

However these networks and the other related models, in their standard 
formulation, are unable to describe the behaviour of a material subjected to a 
creep-test, which consists in keeping a sample at a constant stress till it fails. 
Creep-tests are widely used by engineers in order to estimate the sample life 
time as a function of the applied stress. A modified fuse network, which takes 
into account the Joule effect in the fuses, has been proposed to explain the 
finite life time of a sample subjected to a constant stress^]. However also this 
model does not explain the recent experimental results on micro crystals HIT 



gels[]T2| and heterogeneous materials |13[. These experiments show that the 
life time r of a sample, subjected to an imposed stress P is well predicted by 
the equation 

/ pdy(d-l) \ 
r = r Q exp (a — p(M _ 2) j (1) 

where r D is a constant, T the surface energy, Y the Young modulus, k 
the Boltzmann constant, a a constant which depends on the geometry, T e ff 
is an effective temperature and d the dimensionality of the system. The 
main physical hypothesis behind eq.p] is that the macroscopic failure of a 
material is produced by a thermal activation of micro cracks|n, [T5|]. In the 



original Pomeau's theory |14| T e ff of eq.|TJ coincides with the thermodynamic 



temperature T while experimentally T e ff » TflTJJ, 13]. To explain this 



results it has been supposed that the disorder of the material "amplifies" the 
thermal noise JITJ . 



The purpose of this letter is to show that adding to the spring network 
a noise, which plays the role of a temperature, it is possible to reproduce 
the behaviour of a material subjected to a creep-test. Specifically the func- 
tional dependence on the applied stress of the network life time is consistent 
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with that observed in recent experiments, that is with eq.|I|. Furthermore we 
can prove that, for such a process, the noise is "amplified" by the network 
disorder. 

In order to study the influence of noise on macroscopic failure we have 
chosen the simplest spring network, that is the so called democratic fiber 



bundle model, proposed a long time ago by Pierce [|T(J to study cable failure. 
This model, widely studied in the quasi static regime |], ||, || || , is equivalent 
to N springs in parallel subjected to a total traction force F. Specifically we 
have numerically studied the model by using the following rules: 

Al The external applied force F produces a local stress fi on each fiber. F 
is democratically and completely distributed in the net: F = Yui=i fi 

A2 The local stress /, on the i-th fiber produces a local deformation ej. 
Being the media elastic, stress and deformation are linked by Hook's 
law: 

fi = Y ■ (2) 

where Y is the Young modulus, which is assumed to be the same for 
all the fibers : Y t = Y. 

(c) 

A3 The strength of each fiber is characterized by a critical stress / 4 : if at 
time t on the i-th fiber local stress fi(t) is greater than critical stress / 4 
the fiber cracks, and his local stress falls to zero at time t + 1. Further, 
we assume that in this process some energy is released proportionally 
to the square of local stress fi, for sake of simplicity we will assume 
6j = ff. Critical stress f^ is a realization of a random variable that 
follows a normal distribution of mean f^ and variance KT^. 

= /M + N d {KT d ) (3) 

We call Nd disorder noise. 

A4 Each fiber is subjected to an addictive time dependent random stress 
Afi(t) which follows a zero mean normal distribution of variance KT: 

Afi(t)=N T (t,KT) (4) 

being t the time. We call Nt thermal noise. We assume that A/j(t) 
is a white random process,which is independent in each fiber, i.e. the 
correlation function E [A/j(£i) ■ A/ ? (t 2 )] = if t 1 ^ t 2 or i ^ j. 
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The first three items are those used in the standard formulation of fiber 
bundle model. A new feature, which is similar to a thermal activation process, 
is introduced in [A4] to explain the dependence of failure time on a constant 
applied stress. 

The model has the following properties. We see from [Al] that there 
exists a long-range interaction between the fibers: indeed if a number n(t) of 
fibers are broken at time t, the local stress on each of the remaining fibers 
will be: 

that is, under the specified boundary conditions, the breaking of some fibers 
produces an increase of local stress on the other ones. Item [A3] models the 
heterogeneity of media. The assumption that all the disorder in the model 
appears in the strength distribution rather than in the elastic constants may 
be argued by noticing that the effective elastic constant of a single fiber is 
essentially the average of the local elastic constant along the fiber, while 
the strength is determined by its weakest point 0. If KT = the model 
is reduced to the standard one: in this case the applied force is increased 
linearly F = A p t from zero to the critical value F c needed to break the whole 
network. The sample life time r is equal to F c /A p for any value of the slope 
A p . 

It is quite clear that when KT = and a constant force F is applied the 
system breaks in a single avalanche only if F is large enough otherwise it will 
never break. For example if KT^ = all the fibers are strictly equal and if 
there is no thermal activation ( i.e. KT = 0) the system breaks in a single 
step when F = Nf( c >. 

If KT ytz o then the system can break also at constant imposed stress. We 
have numerically studied the behaviour of the model as a function of F, KT 
and KT&. In the following we assume N = 1000 and / (c) = 1. For each set 
of parameter values we have repeated the numerical simulation at least ten 
times to estimate the scattering of the results for different realization of the 
noise. This corresponds to the scattering of symbols in the figures. We call 
event the simultaneous breaking of several fibers and event size the number 
s(t) of fibers which crack. The energy e associated to an event is the sum of 
the energies released by the fibers which crack, that is : 
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The cumulated energy E(t) is the sum of e from to t. When KT ^ and F 
is constant the event statistics is quite similar to that observed at KT = for 
F increasing linearly in time Jl], 0, |5|] . As an example of the system response 
to a constant force F = 540 with KT d = 0.005 and KT = 0.007, we plot 
in fig. la the distribution N(e) of e and in fig. lb E(t) as a function of the 
reduced parameter (r — t)/r (notice that because of the constancy of F the 



only control parameter is time). In agreement with experiments fT9|, [13 
we find that e is power law distributed and that the cumulated energy E 
has a power law dependence on (r — t)/r for t — > r. More details on this 
statistical features of the model for different driving forces will be given in a 
longer report. Here we want to focus on failure time r of the network in the 
case of a creep test (constant F). 

We first keep KTd = and and we study the evolution towards failure 
for various KT and F. The results are summarized in fig. 2. When KT ^ 0, 
we observe (see fig. 2a) that failure time r as a function of follows an 
exponential law for any fixed value of KT^. 



exp 



F 



(7) 



where Fq is a fitting parameter. Further, looking at fig. 2b we notice that the 
failure time r depends on thermal noise KT as follows: 

r ~ exp {4f) (8) 

where A is a fitting parameter. 

We notice that these results are similar to the prediction of Pomeau's 
theory, that is eq.l with d = 2. Because of this analogy with eq.[TJ, we are 
now interested in studying the dependence of failure time r on disorder noise 
(i.e. on KTd): to this aim, we have done simulations keeping the thermal 
noise variance KT fixed to constant (not zero) values. Looking at fig. 3, we 
observe the following facts: 

Bl failure time r decreases following a power law as the disorder noise 
variance KTd increases, that is the more the media is heterogeneous, 
the smaller is failure time (see fig. 3a). 

B2 as KTd increases, the absolute difference between failure times r for 
different values of thermal noise variance KT decreases, that is failure 
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time r becomes less sensitive to the effective value of thermal noise, as 
shown in fig. 3b. 



Thus one may conclude that disorder noise amplifies the effect of ther- 
mal noise and reduces the dependence of r on the temperature; these results 
allow us to explain the recent experimental observations. Let us recall that 
experiments on microcrystals |TT|, gels [0 and macroscopic composite ma- 



terials reveal that dependence on P of the sample life time is very well 
fitted by eq. ([!]). However calculations [13] show that thermal fluctuations are 
too little to activate the nucleation of microcraks in the times r measured 
in the experiments. It had been measured that the temperature needed to 
have the measured life-times r should be of the order of several thousands 
of Kelvin (3000if for wood [|I1J ). Experiments show that the life time r 



of very heterogeneous materials [13], is (in the limit of experimental errors) 
independent of T while the life time r of quite homogeneous materials as 
microcristals |Tl]] strongly depends on T, as predicted by eq. ([[]). To explain 
these results has been supposed that disorder amplifies the thermal noise so 
that the nucleation time of defects becomes of the order of the measured 
ones and that the lifetime r of the sample depend on the heterogeneity of 
the media. These hypothesis are now well verified by the numerical results 
Bl and B2 of the fiber bundle model. 

As a conclusions we have shown that adding a white noise, which plays 
the role of temperature, to the fiber bundle model we can reproduce the 
recent experimental observations on the dependence of sample life time on the 
applied stress. Furthermore using such a simple model we can prove that the 
disorder of the system amplifies the thermal noise in crack nucleation process. 
This explains quite well why the experimentally measured T e jf is more close 
to the thermodynamic temperature in microcrystals than in heterogeneous 
materials. Similar conclusions about a disordered induced high temperature 
in nucleation processes have been reached in other disordered systems such 
as foams |Tj| . This is quite interesting because it seems to be a quite general 
properties of disordered systems where a thermal activated processes with 
long range correlation may be present. 

We acknowledge useful discussion with Y. Pomeau. One of us (R.S.) 
thanks "Le Laboratoire de Physique de l'E.N.S.L." for the very kind hospi- 
tality during his visit in Lyon. 
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Figure 1: (a) Distribution of energy N(e) symbols are numerical measures, 
while continuous line represents a power law fits N(e) ~ of small events 
with e < 10. (b) E(t) as a function of the critical parameter — ; circles are 
numerical measures, while the line represent the best fit, near t ~ r,with the 
function E(t) ~ (— ) 7 . The parameters used are KT = 7 10~ 3 , KTd = 
5 lO" 3 and F = 540. 
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1/KT 

Figure 2: Failure time r of an homogeneous net (KT^ = 0) in a creep test, 
(a) r as a function of the normalized force / = F/N for several values of 
thermal noise variance KT [(o) KT = 9 1(T 3 , (<})KT = 12 1(T 3 , (*)KT = 
15 10~ 3 ]. (b) r as a function of KT for several values of f [ (*)/ = 0.63, (0)/ = 
0.6, (A)/ = 0.56, (o)/ = 0.53]. Continuous lines in (a) and (b) are the fits 
with eq.7 and eq.8 respectively. 
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Figure 3: Failure time r of an heterogeneous network (KT d ^ 0) in 
creep test, (a) r is plotted as a function of normalized force / = F/N, 
at KT = 6 10" 3 , for different values of KT d [ (A)KT d = 9.8 1(T 3 , (o) KT d = 
8 1(T 3 , (0)KT d = 7 KT 3 , (*)KT d = 6 KT 3 ]; continuous lines represent the 
best fits with eq. 7. (b) function of disorder noise variance KT d at 

F = 540 for several values of KT [KT X = 7 10" 3 , KT 2 = 8 KT 3 , KT 3 = 1(T 2 ]. 
Dashed lines represent the data fits . 
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